Controlling phase separation of a two-component Bose-Einstein condensate by 

confinement 
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We study the effect of kinetic energy on the phase separation and phase transition of a two- 
component Bose-Einstein condensate in the presence of external confinement. The commonly ac- 
cepted condition for the phase separation, giig22 < gi2 where gn, (722, and 1712 are the intra- and 
inter- component interaction strengths respectively, is only valid when kinetic energy is negligible 
and external confinement is nonexistent. Taking a d-dimensional infinitely deep square well po- 
tential of width L as an example, a simple scaling analysis shows that regardless of the condition 
911322 < 312, if d = f (d = 3), phase separation will be suppressed as L — >■ (L — > 00) and if d = 2, 
the width L is irrelevant but again phase separation can be partially, or even completely suppressed. 
Moreover, the miscibility-immiscibility transition is turned from a first-order one into a second-order 
one when kinetic energy is considered. All these results carry over to d-dimensional harmonic poten- 
tials, where the harmonic oscillator length S^ho plays the role of L. Our finding provides a scenario 
of controlling the miscibility-immiscibility transition of a two-component condensate by changing 
the confinement, instead of the conventional approach of changing the values of the g's. 

PACS numbers: 03.75.Mn, 05.30.Jp 



I. INTRODUCTION 

Phase separation is a ubiquitous phenomenon in na- 
ture [J, Q • A most prominent example familiar to ev- 
eryone is that oil and water do not mix. Besides that, 
the phenomenon of water in coexistence with its vapor 
can also be understood as a type of phase separation [31 ■ 
In general, two phases mix or not depending on which 
configuration minimizes the energy or free energy of the 
whole system. With the realization of Bose-Einstein con- 
densation in ultracold atomic gases, another example 
of phase separation is offered by two-component Bose- 
Einstein condensates (BECs) [1-11]. In such a system, 
phase mixing or separation means the two condensates 
overlap or not spatially, which correspond to different in- 
teraction energies. A widely accepted condition for phase 
separation, which is based on the consideration of mini- 
mizing the interaction energy [lol | , is given by 

(1) 



512 > V511522- 

Here gn and 522 are the intra-component interaction 
strengths of components 1 and 2, respectively, while gi2 is 
the interaction strength between them [ll|. This condi- 
tion is intuitively reasonable since if the inter-component 
interaction is too strong, the two components would 
like to get separated from each other. Experimentally, 
controlled miscibility-immiscibility transition of a two- 
component BEG based on the idea of adjusting the values 
of the g's using Feshbach resonance and so as to ge t 
satisfied or not has been demonstrated recently |12l \lc 
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Now the point is that though the condition above is 
very appealing in its simplicity and usefulness, it has 
great limitations. In its derivation, the condensates are 
assumed to be uniform and the kinetic energy associ- 
ated with the boundary /interface layers is neglected. The 
problem is then reduced to minimizing the total interac- 
tion energy, or more specifically, to weighing the inter- 
component interaction against the intra-component in- 
teraction. This approximation is legitimate if the widths 
of the boundary/interface layers are much smaller than 
the extension of the condensates, or in other words, if 
the boundary/interface layers are well defined. How- 
ever, this condition is not necessarily satisfied in all 
circumstances. Actually, some simple scaling analysis 
may tell us when it will fail. Consider a condensate 
trapped in a d-dimensional container of size L. The 
characteristic (average) density of the condensate is on 
the order of L~'^. According to the mean-field (Gross- 
Pitaevskii) theory, the healing length of the condensate, 
which determines the widths of the boundary/interface 
layers, will be on the order of L'^^^ Thus we 

see that in one and three dimensional cases, it makes 
sense to say boundary/interface layers only in the lim- 
its of L — >■ 00 and L — >■ 0, respectively. In the opposite 
limits, the "boundary/interface" layers overtake the con- 
densates themselves in size, which signals that the kinetic 
energy will dominate the interaction energy and should 
no longer be neglected. The two dimensional case is more 
subtle in that the widths of the boundary/interface layers 
scale in the same way with the sizes of the condensates, 
which at least means that the kinetic energy should not 
be neglected a priori. 

The analysis above indicates that the kinetic energy is 
likely to play a vital role in determining the configura- 
tion of a two-component BEG. Moreover, we note that 
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the kinetic energy acts against the inter-component in- 
teraction. The latter is responsible for phase separation 
while the former tries to expand the condensates and 
thus favors phase mixing. Therefore, it is expected that 
phase separation can be suppressed by the kinetic en- 
ergy in some circumstances even if the condition ^ is 
satisfied Ts']. Notably, according to the argument above, 
the significance of the kinetic energy can be controlled 
by changing the size of the container. That is, the phase 
mixing-demixing transition can be controlled by a geo- 
metrical method, instead of the mechanical method of 
changing the values of the q's, which is based on ([1]) and 
is demonstrated in Refs. |12l. [l3|. 



II. A TWO-COMPONENT BEC IN AN 
INFINITELY DEEP SQUARE WELL POTENTIAL 

The considerations above have led us to investigate 
the scenario of suppressing phase separation in a two- 
component BEC by kinetic energy. We will start from 
the simplest and most generic case of a two-component 
BEC in a d-dimensional infinitely deep square well po- 
tential (of width L). The Dirichlet boundary condition 
implies that the condensate wave functions must be non- 
uniform and the kinetic energy is at least on the order 
of L~^. On the contrary, inside the well, the potential 
energy is zero. Therefore, we have a pure competition 
between the kinetic energy and the inter-component in- 
teraction energy, if the intra-component interactions are 
set zero [note that in this case, condition ((TJ is satisfied] . 
In this simplest model, in all dimensions {d = 1, 2, 3), we 
do observe that phase separation can be completely sup- 
pressed by the kinetic energy in some regime. Of course, 
different dimensions have different features. But all these 
effects and features carry over to the more realistic case 
of d-dimensional harmonic potentials. 

In the mean-field theory and at zero-temperature, 
the energy functional of a two-component BEC in a d- 
dimensional infinitely deep square well potential — 
[-L/2, +L/2Y is of the form 

E[4>,,^,]^ [ dr[ J2 ^IV^.P 

+ 1 E 9a^N^Np\i;^\^\i;^A. (2) 

a, 13=1,2 ■' 

Here the two condensate wave functions are normalized 
to unity rf'^|'0i,2p = 1, and ^pi^2 = on the boundary. 
Note that throughout this paper we are only concerned 
with the ground configuration of the system, therefore 
all the wave functions can be taken to be real and posi- 
tive. The parameters gu, 522, and 512 = 521 are the ef- 
fective intra- and inter-component interaction strengths. 
Finally, and mi^2 are the atom numbers and atom 
masses of the two species, respectively. Now we should 
note that for an arbitrary set of parameters, in the ground 
configuration, almost definitely, the two wave functions 



do overlap but do not coincide with each other (this can 
be easily understood in terms of the Gross-Pitaevskii 
equations for "01. 2)- In this case, it is far from trivial 
to distinguish phase separation and phase mixing. A 
method proposed in ^IJJ is to consider the centers of mass 
of the two condensates: 



dr\^pc 
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1,2. 
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This idea is motivated by the observation that in some 
regime, both the two condensates are symmetric with re- 
spect to the origin while in other regime, both of them are 
asymmetric with respect to the origin, and more impor- 
tantly, they are shifted in opposite directions fl4l|. Ap- 
parently, the former case is with r„ji ~ r„j2 = and it is 
appropriate to call it phase-mixed while the latter case is 
with r,„i 7^ ^ 7Vi2 and it is appropriate to call it phase- 
separated. Therefore, the offset between the two centers 
of mass fmi — r'm2 can serve as an order parameter for 
the miscibility-immiscibility transition of the system. 

Though this order parameter works well for a general 
case, we will not use it much in this paper. Actually, in- 
stead of studying a general case, we shall focus on the 
symmetric energy functional case, i.e., the case when 
mi = 1712 = m, Ni = N2 = iV, and gu = g22- The 
reason is that this special case not only captures all the 
essential physics, but also has an extra merit. That is, 
now it is possible to have ipi = "02, which corresponds 
to a completely mixed configuration. Therefore, in this 
special case, an appropriate order parameter is the over- 
lap between the two condensate wave functions (or more 
precisely, 1 — 77, if phase separation is concerned): 



T] ^ df?/;ii/'2, 
Jn 



(4) 



which takes values between and 1. If 77 ^ 1, it would be 
fair to say the system shows phase separation. Otherwise, 
if 77 is close to 1, or more precisely if 1 — 77 ^ 1, it would 
be fair to say the system shows phase mixing. In the 
intermediate case, the system is partially phase-separated 
and partially phase-mixed. 

Now make the transform V'i,2(^) — i~''^^0i,2(^) with 
f = Lx. Then J^^dx\(f)i^2\^ — 1 and 0i_2 = on the 

boundary of Oq, where VIq = [—1/2, +1/2]''. In terms of 
the rescaled wave functions 0i,2, = /q^^ dx(j)i(t>2, and 
the energy functional ([2]), under the assumption above, 
can be rewritten as 
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with the reduced dimensionless parameters /3ij defined as 

r, «,i = l,2. (6) 
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FIG. 1. (Color online) (a) The overlap factor 77 as a function of 
the reduced parameter /3i2 [see Eq. ^] in different dimensions 
(infinitely deep square well potential case, gn — g22 ~ 0). 
Note that for all values of d, there exists a critical value /?J2 7^ 
0, below which 77 attains its maximal possible value 1. (b) a 
schematic plot of rj versus the width of the square well in 
different dimensions. Note the counter-intuitive fact that in 
the three dimensional case (d = 3), the stronger we squeeze 
the system (the smaller L is), the stronger phase separation 
is (the smaller rj is). 

These parameters are measures of the importance of the 
interactions. In the curl bracket, the coefficients of the 
kinetic terms are constant, yet the coefficients of the in- 
teraction terms (the /3's) scale with L as L^^'^. This fact 
has some important consequences. If d = 1, there are 
two different limits. In the limit of L — > 00 (loose confine- 
ment), the kinetic terms are dominated by the interaction 
terms and thus the ground state can be determined by 
simply minimizing the interaction energy. In this limit, 
the textbook analysis is valid and we have phase separa- 
tion if condition ((U is satisfied or phase mixing otherwise. 
In the opposite limit of L — >■ (tight confinement), the 
kinetic terms will dominate and the two rescaled wave 
functions can be well approximated by the ground state 
of the square well potential, i.e., 4'i,2ix) — cos(7ra;). In 
this limit, phase separation will be suppressed whatever 
the values of the g's are, even if ([T]) is fulfilled. The three 
dimensional case is the inverse of the one dimensional 
case. In the limit of L — ^ 0, the kinetic terms are neg- 
ligible and the criterion of phase separation ([1]) is valid. 
In the other limit of L — >■ cx) , the kinetic terms dominate 
and phase separation is suppressed regardless of the con- 
dition ([T]). The two dimensional case is another story. 
The parameter L simply drops out in the curl bracket. 
It is no use to adjust the width of the well to enhance 
the importance of the kinetic energy or the interaction 
energy relatively. The kinetic and interaction energies 
should be treated on an equal footing, which means the 
analysis leading to criterion ([T]) may be invalid. 

We have checked all these predictions numerically. 
Note that on the problem of phase separation, the intra- 
component interactions are on the same side as the ki- 
netic energy — they both try to delocalize the conden- 
sates. Therefore, to highlight the effect of kinetic energy, 
we shall set gn = 522 = (/3ii — P22 = 0), so that the 
kinetic energy is the only element acting against phase 



separation. As we shall see below, this special case also 
admits a simple analytical analysis. 

We have solved the ground state of the system in all di- 
mensions for a given value of /3i2 fl6} . The overlap factor 
77 is plotted versus /3i2 in Fig. [T^. We observe that in all 
dimensions, there exists a critical value of /3i2 (denoted 
as /3f 2 ) : below which the two condensates wave functions 
are equal {-q = I). That is, for (3i2 < /JJjj phase separa- 
tion is completely suppressed. Above the critical value, 
phase separation develops (77 < I) as /3i2 increases, but is 
still greatly suppressed for a wide range of value of /3i2. 
It should be stressed that though in Fig. the curves 
of 7/ — /3i2 are qualitatively similar to each another for all 
values of d (the plateau of 77 = 1 is always located in the 
direction of (3i2 0), the curves of 77 — L will be quite 
different. The reason is that /3i2 oc i^^'^. Figure [Dd is 
a schematic plot of 77 versus L in all the three cases. It 
shows that 7; as a function of L is monotonically decreas- 
ing, constant, and monotonically increasing in one, two, 
and three dimensions, respectively. This means that to 
suppress phase separation, in one dimension we should 
tighten the confinement, in three dimensions we should 
loosen the confinement, while in two dimensions it is use- 
less to change the confinement. Overall, Fig. [T] confirms 
the initial conjecture that kinetic energy can suppress 
phase separation. 

As a hindsight, we can actually understand why phase 
separation can be suppressed in the limits of i — > in 
one dimension and i — >■ 00 in three dimensions. Consider 
two different configurations. The first one is a phase- 
separated one — the two condensates occupy the left and 
right halves of the container separately. The second one 
is a phase-mixed one — the two condensates both occupy 
the whole space available and thus overlap significantly. 
Compared with the first configuration, the second one 
costs more inter-component interaction energy which is 
on the order of L^"^, but saves more kinetic energy which 
is on the order of L~^. The second configuration (phase- 
mixed) is more economical in energy in the limit of L — 
and L — > 00, in the cases of d = I and d = 3, respectively. 
The case of d = 2 is more subtle and which configuration 
wins depends on parameters other than L. 

A remarkable fact revealed in Fig. [T] but not so obvi- 
ous in our arguments is that in the symmetric case with 
= /?22 = 0, 77 = 1 for /3i2 < /3i2i which is on the order 
of unity. This is a stronger fact than 77 — >■ 1 as /3i2 — > 
as we argued. Actually, the general observation is that 
for Pii — (322 > 0, ?7 = I for /?i2 smaller than its critical 
value I3i2i which is larger than /3ii. This fact has rich 
meanings. On the one hand, it demonstrates that the 
kinetic energy is very effective — phase separation can be 
completely suppressed by it even if j3i2 > /3ii = 1322, i-C, 
when ID) is satisfied. On the other hand, it strongly in- 
dicates that as /3i2 crosses the critical value, the system 
undergoes a second order phase transition which can fit 
in the Landau formalism. The picture is that the ex- 
change symmetry 0i -f-)- 02 of the energy functional ([5]) 
is preserved for /?i2 < /3x2i but is spontaneously broken 
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as /3i2 surpasses /3J2- 

We have been able to prove the first point rigorously 
on the mathematical level (see Appendix R)) . However, 
it is also desirable to develop a physical understanding 
of the two points. This can be achieved by studying 
a two-component BEC in a double-well potential (see 
Appendix or using a variational approach [13] . We 
note that in the limit of P12 — >■ 0, (/)i,2 both converge to 
the (non-degenerate) ground state of a single particle in 
the [—1/2,-1-1/2]'' infinitely deep square well. As /3i2 is 
turned on, the two wave functions are deformed and ex- 
cited states mix in. Because the energies of the excited 
states grow up quadratically, we cutoff at the first excited 
level and take the following ansatz for the two condensate 
wave functions 



CQipo+Ciifi, (f>2 = CQipo - Clipi. 



(7) 



Here ipo is the ground state, while (pi is one of the possibly 
degenerate first excited states. The coefficients co,i are 
real and satisfy the normalization condition Cq + = 1. 
Obviously, complete phase mixing would correspond to 
ci = while partial phase separation to ci ^ 0. Our nu- 
merical simulations indicate that (this is also supported 
by the variational approach itself, see Appendix [C]) in the 
two dimensional case, when phase separation occurs, the 
two condensates are shifted either along x or y direction; 
in the three dimensional case, when phase separation oc- 
curs, the two condensates are shifted either along x or y 
or z direction. This fact motivates us to choose (fii in the 
following form 

d — 1 : ipi — wi{x); (8a) 
d = 2 : (^1 = WQ{x)wi{y) or wi{x)wo{y); (8b) 
d = 3 : (^1 = wo{x)wo{y)wi{z) or wo{x)wi{y)wo{z) 

or wi {x)wo {y)wo (z) , (8c) 

where wo{x) ~ -\/2 cos(7ra;) and wi(a;) = \/2sin(27rx) are 
the ground and first excited states of a single particle in 
the one dimensional [—1/2,-1-1/2] infinitely deep square 
well potential. Substituting Eqs. d?]) and ([S]) into (O, we 
get the reduced energy functional E — E j (NT? jml?^ as 

d=\: E{ci) ^ {3n^ - 6^12)0? -|- b/Sucf + const; 

d = 2 : E{ci) = (^3tt^ - y/3i2^ cj + ^-Pv2c\ + const; 

~ / 45 \ 45 

d = 3: E{ci) = ( Btt^ - — /?i2 j + ^^2^ + const. 

These are nothing but the Landau's expression of the free 
energy in a second-order phase transition, with ci play- 
ing the role of the order parameter here. We immediately 
determine the critical values of /3i2 by putting the coeffi- 
cients of c{ to zero. Specifically, P12 = and ^ 
for d = 1, d = 2, and d = 3, respectively. These values 
agree with those extracted from Fig. [T] very well. The 
relative errors are within 1%, 9%, and 19%, respectively. 



The deviation increases with d because in higher dimen- 
sions, the degeneracy of the excited states increases and 
the two-mode approximation in Q becomes less accu- 
rate. In the expressions of E, we can actually see how 
the kinetic energy suppresses phase separation. The term 
37r^cf comes from the kinetic energy difference of the two 
modes </9i,2- Without this term, the critical value /3f2 
would be zero instead of being finite. 

For a general case without the exchange symmetry 
(j)i -f-)- (j)2^ the appropriate order parameter is no longer 
rj but rvni — fm2- However, the second order transition 
picture still holds. Specifically, f,„i = = f,„2 for /3i2 
smaller than some critical value /3J2 which is larger than 
■\//3ii/322- Overall, this asymmetric case is more involved 
than the symmetric case above because there are more 
parameters. Hopefully, a systematic study will be pre- 
sented in a follow-up work. 



III. A TWO-COMPONENT BEC IN A 
HARMONIC POTENTIAL 



So far, we have focused on the ideal case of infinitely 
deep square wells. Experimentally, it is harmonic po- 
tentials that are most readily realized. Therefore, it is 
necessary to see whether analogous results hold for har- 
monic potentials. One concern is that the extra potential 
energy may blur the picture. However, after some similar 
rescaling, we shall see that all the results persist. 

The energy functional of a two-component BEC in a 
d-dimensional isotropic harmonic potential is 

" ^ a=1.2 a=\.2 

+ Y (fflll^ll' + 522|V'2|' + 2gi2|Vl|'|V'2|') |(.10) 

Here again we have assumed equal mass and equal num- 
ber for the two species. The two condensate wave func- 
tions are normalized to unity, i.e., /dr|V'i,2|^ = 1- Now 
make the transform '^\.2(j) = ^^o^"^ 4>\,2(x) with r = 
£,hoX, where £,ho — \/nJmiUd is the characteristic length 
of the harmonic potential. We have then J dx\(j)i^2\'^ — 1. 
In terms of (/>i,2, the energy functional can be rewritten 
as 



E 



+ 2 (Ail'^i 



-/322|</'2|^+2/3i2|(/.l|'|(/.2 



,(11) 



Here the reduced interaction strengths are defined as 



ft: 



ft2 



{d-2)/2 



1,2. 



(12) 



We now have a similar situation as before. The im- 
portance of the interactions can be changed by chang- 
ing the value of ^ho, which plays the role of L in our 
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FIG. 2. (Color online) (a)-(c) The overlap factor r; as a func- 
tion of the reduced parameter /3i2 [see Eq. (|12|) ] in different 
dimensions (isotropic harmonic potential case, gn = (?22 = 
(712/1.05). Note that for all value of d, there exists a criti- 
cal value /?i2 7^ 0, below which 77 attains its maximal possible 
value 1. (d) a schematic plot of 77 versus the characteristic fre- 
quency LOd of the harmonic potential in different dimensions. 



previous example. The interactions will be negligible 
if d = 1 and i^ho — or if d = 3 and ^/lo 00. 
In this case, the rescaled wave functions (/>i^2 will be 
close to the ground state of the harmonic oscillator, i.e., 
(/)i_2 — TT"'^/^ exp(— and phase separation is sup- 
pressed regardless of the values of the t/'s. The interac- 
tions will become significant if d = 1 and S^ho — >■ c» or 
d = 3 and ^ho —^0. In this case, the kinetic energy can 
be neglected and we enter the Thomas-Fermi regime. In 
this regime, the criterion (jlj will be a faithful one for 
phase separation. 

We have verified these predictions numerically. In 
Fig. [21 we have shown the overlap factor r] = J dx<j)i(t>2 
versus the reduced inter-component interaction strength 
/3i2 in all dimensions (with gn = (722 = 312/I.O5). Again, 
we see that phase separation is completely suppressed for 
/3i2 below some critical value I3i2- 

Let us now consider the possibility of experimentally 
observing the immiscibility-miscibility transition by ad- 
justing the confinement, e.g., the frequency ujd- In cold 
atom experiments, the harmonic potential is often of the 
form V{f) = ^m[uj\{x'^ +y'^) +^lz^]- To get a three di- 
mensional isotropic potential, we set ujj_ — uj^ ■ An effec- 
tively one (two) dimensional potential can be obtained in 
the limit of uj± 3> oJz (w± <^ ojz). For these three differ- 
ent geometries of the potential, the interaction strengths 
(the g's) relate to the s-wave scattering lengths (the a's) 
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LOd^^z ^ w^, d = 3; 



(13a) 



^, Ud=uj^^uz, d=2; (13b) 

(13c) 



gij = 2haijLL)±, uJd = ujz d — 1 



Using Eqs. (|T2t and (|T3)) . we can study the possibility 
of tuning /3i2 across the critical value f3i2- We study 
each case individually (the mass m is taken to be that of 
23Na): 

(i) d = 3. Suppose N — lO'*, ai2 — 40 ob. The critical 
value of LJd is 27rx 560 Hz, which can be covered in current 
experiments. 

(n) d — 2. Suppose N = 10^, ai2 — 40 ae, and the 
transverse frequency uj± = 27r x 2.6 Hz. The critical value 
of the longitudinal frequency uj^ is 27r x 140 Hz, which is 
realizable in current experiments [l2j . 

(iii) d = 1. Suppose N = 2 x 10^, ai2 = 40 ae, and 
the transverse frequency uj± = 27r x 130 Hz. The critical 
value of the longitudinal frequency is 27rx 19 Hz, which 
is realizable in current experiments. 

Here the number of atoms is one or two orders smaller 
than its typical value in experiments. This explains why 
the criterion ([1]) is a reliable one in the experiments in 
[13, [3 • They work in a regime where the kinetic energy 
is indeed negligible. However, with the advance of imag- 
ing techniques, hopefully future experiments can work 
with a relatively small number of atoms and observe the 
miscibility-immiscibility transition by changing the con- 
finement. 



IV. CONCLUSIONS 

To conclude, we have demonstrated that kinetic energy 
can play a vital role in determining the configuration of a 
two-component BEC. It renders the empirical condition 
of phase separation giig22 < <?i2 insufficient and it also 
modifies the picture of phase separation. To be specific, 
phase separation can be completely suppressed even if 
this condition is fulfilled. Moreover, the phase mixing to 
phase separation transition is now known to be a second- 
order, continuous transition instead of a first-order, dis- 
continuous one as in the usual view. From the experimen- 
tal point of view, our results may provide a new scenario 
of controlling the transition of phase mixing- demixing of 
a two-component BEC. Instead of adjusting the interac- 
tion strengths, one can just change the confinement, the 
characteristic size of the container. 
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Appendix A: Rigorous justification 

Here, we consider the energy functional as 

|4 I o IJ, |4 , na \A |2 



:|V02|' 



(Al) 



1 



+ 2 (/^"I'^ll" + /^22|02|' + 2/3i2|0in02|' 



where = [-5,5]'' (d = 1,2,3), = ^22 = Let 
(pg be the unique positive ground state of the energy 
functional Es[(j)] = E[(j),(j)], and be the corresponding 
chemical potential. The functions 0i,2 are normalized to 
unity by the usual L^-norm. Let be the pos- 

itive ground state of (jAip . For /3i2 < /3, E[y/pi, y/p2] 
(£1 = l^ip, p2 = |(?!'2p) is strictly convex in (pi,/92) 
[ill, [1^, and the positive ground state is unique, i.e., 
4'i = 4'2 — 't'g^ V = 1- We are going to prove that there 
exists a critical value f3i2 > P such that when /3i2 < /?i2: 
there holds <jT{ = 0f — 4'gi rj — \. From now on, 
we concentrate on the case of /3i2 > /3 and assume that 
P12 = /3 + /3',0< /?'<!. Simple calculation shows that 



E\ 



dx 



+ {P + Pl2)\4>a?\<t>l.-<t>g? + 



a=l,2 

/3-/3: 
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2 VITil ■ ITZI -iV^Sl 

Making use of the Euler-Lagrange equation of (f)g , 



(A2) 



denoting = 



, — (pg (ck = 1,2), and noticing 
ill Sail 2' obtain 



E[ct>l<l,l]-E[^g,(^g]^ j dx^/l|Ve„| 



(/3 + /3i2)|0,ne„p-^(|0?:i^-|0, 



2fig(j)gea 



10: 



IP 



Now, the operator Lg = — ^V^ + {f3 + Pi2)\4'g\'^ admits 
eigenvalues as /ig < /ii < /i2 < • • • , and the eigenfunc- 
tion (pg corresponds to /ig, Wk G Hq with ||wfc||2 = 1 
corresponds to fik (k > 1) . The reason 0g is the ground 
state comes from the positivity of 4>g and the uniqueness 



of the positive ground state of Lg. Expand as = 
c'^'l^g + EZi^>k, then (c^)^ + Er=i K\' = 
Cg — J^^ BaCpg — —5 II Call 2 we Can derive that 

|i|Ve„P + {/3 + Pl2Mg\^K\^ + ^flgCPgC^^ 

00 



k=l 



>(/ii-M.)(l|e.||^-(c^)') 
-(Mi-Mg)l|ea||^(l-||<||^/4). 

Now, firstly, we need a lower bound for fii — ^ig, the so- 
called fundamental gap, which has been solved recently 
by A. Ben and C. Julie 2,0]. Using equation (IA2p . apply- 
ing elliptic theory with convex domain fio, it is easy to 
verify that (pg € H'^iyio) and hence belongs to C^'''{Q,q) 
(0 < 7 < ^) by Sobolev embedding. Approximating fig 
by convex domain $7^ (with smooth boundary) and ap- 
plying Schauder estimates, we shall have (pg S C'^''^{^^) 
and there exists some c > such that |(/)gp -|-c|a;p is con- 
vex (as Hessian matrix of \4>g\'^ is bounded by Schauder 
estimates). Hence, we can apply the results in Ref. ^2^ 
to get (Z?j is the diameter of Jig) 



Mi - m; 



> 



37r^ 



(A3) 



where /i^ and /if are the first and second eigenvalues, 
respectively, of Lg in 7fg(r2e). By Min-max principles, 
letting £ — >■ 0, we have /ig — >■ /ig and /if — /xi. Hence we 
find 



Ml - Mg 



> 



37r2 
D2 



(A4) 



where D is the diameter of ilg [or if we assume fio is 
a convex domain with smooth boundaries, (|A4[) follows 
directly] . 

Secondly, we have l|eal|i < dfdf/'^P + Ic^gH = 2. 
Thirdly, we would like to derive L°° bounds of 4>g and 
0^. The Euler-Lagrange equation for 0^ reads as 

M^a^ = -^V^e +/?iei'e +/3i2iei''/'^> (A5) 

with a' 7^ a. For the nonlinear eigenvalues, we have the 
estimates /ig, < 2i;[0f , 0^] < 2£',.[0g], /ig < ^^^.^g] and 
Es [4>g\ can be bounded by choosing any test function (like 
the ground state of —A), which gives Es[(f)g\ < C{1 + f3) 
{C depends on Qq). 

If /3 > 1, considering the point xq G JIq where (pg 
attains its maximum, then A0g(a;o) < and from (jA2l) . 
we have 

fJ-g'Pgixo) >{f3 + l3i2)\(pg{xo)\'^(j)g{xo), 



which gives \\(j)g\ 



< 



< 2C. Similarly, we can 



obtain the L°° bound for 0^ using the Euler-Lagrange 
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equation and < ^ < 4C. Thus, Ug + ^jp^ < 

12C. Combining the three observations above, we get 



1 V2\ 



El 



( Stt^ _ 12/3'C 



Q = l,2 



-0:1121 



which imphes that for < /3' < mini there 
must hold Cq = 0, i.e., rj — 1. 

For /3 g [0, 1], the approach above is not good. In this 
case, we see that /i^ < 4C and fig < 2C. Using Sobolev 
inequahty, in one dimension (d = 1), we can find that 



(A6) 



Similarly, ||0g||^ < V 2C. For two and three dimen- 
sions (d = 2,3), recalling (jA2p and (IA5|) . we can ob- 
tain from elliptic theory and Sobolev inequalities that 
there exist constants Ci , C2 > only depending on 
no such that ||</)^|U < CiUi\\H2 < C2 ■ ll/x^e _ 
mi?4>'L - PM.\H<Lh. and < C2UgCt>g - 

(/3 -I- (3i2)\4>g\^4'g\\2- In two and three dimensions, us- 
ing Sobolev inequality, we have H^^He < C3||V(/)^||2 < 
G-i\f'f^ (C's depends on JIq). Cauchy inequality leads to 

ii0^jioo < c2(m^ + miwi + /3i2ii0^ji6iieiii)> 

and thus 110^11^ < C4 (C4 depends on f2o)- Similarly, 
II '/'s " 



in all dimensions (d = 1,2,3), there exists a constant 
Cq,o depending only on ilp such that ||0q + 0g||L < C'j; 
Similar to the case with /3 > 1, we have 



a=1.2 ^ ■' 



which leads to the conclusion that when f3' < 
min{l, j2g^ }, — 'Pgj i-e., V — ^- In summary, for all 
/3 > 0, if we choose /3f2 = + min{l, ^} > /3, 

then for all < /3i2 < /?i2, we shall have 77 = 1. 



Appendix B: Phase separation as a spontaneous 
symmetry breaking 

Consider a two-component BEC in a symmetric 
double-well potential. Under the two-mode approxima- 
tion, the mean-field energy functional is 

E = -Ja{ll^*al^a2 + ■4^*a2'4^al) - Jf,(V'hlV'b2 + V'mV'&i) 
+ \Ua{\^al\^ + \i^a2\') + ^C/.(|V^w|' + l^bsl") 



+mv'ainv'6ip + i^a2nv'fc2n. 



(Bl) 



Here Ja and Jt, are the hopping amplitudes of the 
two types of atoms, and Ua and C/f, are the intra- 
component onsite interaction strengths, while V is the 



inter-component one. The complex numbers ipai and V'bi 
('0a2 and '062) are the amplitudes of the two condensate 
wave functions on the left (right) trap. They are con- 
strained by the total atom numbers, i.e., |V'aiP + |V'a2p = 
Na and IV'biP + I "0621^ = Nf,. For the sake of simplic- 
ity, in the following we shall assume Ja = Jt = J ^ 0, 
Ua ^ Ub = U > 0, and Na = Nb ^ N. As far as 
the ground state is concerned, it is legitimate to as- 
sume the i/j's real and positive. Therefore, we can write 
ipai = V^ai, '4'bi = V^bi and similarly for other ip^s. 

First assume tunneling is turned off, i.e. J = 0. Let 
Nai ^ + Sa and A^^i = hN - Sb. The energy (HH) is 



E{Sa,6b) = U{6l + Si) - 2VSaSb + const. (B2) 

It is readily determined that ii U > V, the ground state 
is of Sa — Sb — 0. The two condensates are both dis- 
tributed evenly between the two wells, which is a com- 
pletely mixed configuration. If U < V [the counter- 
part of (HI) in the present context], the ground state is 
of {Sa,Sb) — ±(7V/2, iV/2), which corresponds to com- 
plete phase separation — the two condensates occupy the 
two wells separately. Therefore, without tunneling, the 
miscibility-immiscibility transition is a first-order phase 
transition with the critical point being = U. 

Now turn on the tunneling. For the sake of simplicity, 
suppose Sa = Sb ^ S. The energy as a function of the 
order parameter S is 



< C5 (C5 depends on flo). Eventually, we have ^ 




2(U - V)S^ + const 

+ + oiS^) + const. (B3) 

Here we have the familiar Landau formalism for second 
order phase transitions. The coefficient of the quar- 
tic term is positive but the sign of the quadratic term 
changes from positive to negative as V surpasses the criti- 
cal value = U+2J/N. Corresponding, ^ = is turned 
from a minimum to a maximum point and phase separa- 
tion develops. Here we note that the tunneling, the ki- 
netic term in the present context, has two consequences. 
First, the first-order transition is turned into a second- 
order one. Second, the transition point is up shifted from 
U to U + 2J/N. This is reasonable since phase separation 
costs kinetic energy. What presented in Figs. [T] and [5] are 
parallel to these results but in continuum (multi-mode) 
cases. 



Appendix C: Justification of the form of <pi in 
Eq. (HD 

In this Appendix, we show why among all the (degen- 
erate) first excited states, the one in Eq. (jS]) is selected. 
For d = 2, the ansatz more general than Eq. Q is 



CoVO - Cxipx - Cyfy: 



(Cla) 
(Clb) 
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with ip^ = wi{x)wo{y), ^Py = wo{x)wi{y), and cq, c^, 
Cy being some real variables under the constraint Cg + 
+ 4 = 1- Substituting Eq. ([CT]) into Eq. ©, we get 
the reduced energy functional E = E / {N / mL?') as a 
function of Cx,y as 

E[c,,cy] = 27r2 + ^/3i2 + (^Svr^ - ^^^^^ {cl + 4) 

+ f/3i2(c^ + 4)' + ^/3i2C^4- (C2) 

We see that for /3 < /SJj = f""^, the minimum is at 
Cx — Cy = 0. For /3i2 > ^3l2^ the minimum is no longer 
at the origin. However, for a fixed value of + c^, E 
is minimized when the last term in Eq. (|C2|) vanishes or 
when Cx — oi Cy =0. That is why the particular ansatz 
in Eqs. ([7]) and ([S]) is appropriate and enough. We note 
that due to the symmetry of the trap, the reduced energy 
functional is invariant under the transform {cx,Cy) — 
{±Cx,±Cy) and {cx,Cy) {cy.Cx). This symmetry is 
broken when phase separation occurs. 

Similar analysis applies for d ^ 3. In this case, the 



ansatz more general than Eq. (O is 

(t>i = CQipo + Cxfx + Cyipy + c^ipz, (C3a) 
(1)2 = CQipo - Cx'Px - Cytpy ~ Czipz, (C3b) 

with (fx = wi{x)wQ{y)wo{z), ipy = woix)wi{y)wo{z), 
= wo{x)woiy)wi{z), and cq, Cx, Cy, Cz being some 
real variables under the constraint Cq + + + = 1 . 
Substituting Eq. (jC3p into Eq. ([5]), we get the reduced 
energy functional £^ as a function of Cx^y^z as 

E = 3n^~ + f /3i2 + (btt^ - f {cl + 4 + 4) 

+ ^/3i2 {cl + cl+ciy 

+ 1^12 {clcl+ 44 +clcl). (C4) 

We see that for /? < (3^2 — ^""^i the minimum is at 
Cx — Cy = Cz — 0. For /3i2 > the minimum is 
no longer at the origin. However, for a fixed value of 
c^+Cy+cl, E is minimized when the last term in Eq. (|C4p 
vanishes or when two of the three c's are zero. Again, 
we see that the particular ansatz in Eqs. (O and ^ is 
appropriate and enough. 
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